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Introduction. Transverse axisymmetric oscillations of a bimorph with two piezo-active layers, piezoelectric and 
piezomagnetic, are studied. This element can be applied in an energy storage device which is in an alternating magnetic 
field. The work objective is to study the dependence of resonance and antiresonance frequencies, and electromechanical 
coupling factor, on the geometric parameters of the element. 

Materials and Methods. A mathematical model of the piezoelement action 1s a boundary value problem of linear 
magneto-electro-elasticity. The element consists of three layers: two piezo-active layers (PZT-4 and CoFe,Ox,) and a 
centre dead layer made of steel. The finite element method implemented in the ANSYS package is used as a method for 
solving a boundary value problem. 

Results. A finite element model of a piezoelement in the ANSYS package is developed. Problems of determining the 
natural frequencies of resonance and antiresonance are solved. Graphic dependences of these frequencies and the 
electromechanical coupling factor on the device geometrics, the thickness and radius of the piezo-active layers, are 
constructed. 

Discussion and Conclusions. The results obtained can be used under designing the working element of the energy 
storage device due to the action of an alternating magnetic field. The constructed dependences of the eigenfrequencies 
of the resonance and antiresonance on the geometric parameters of the piezoelement provide selecting the sizes of the 
piezo-active layers for a given working frequency with the highest electromechanical coupling factor. 
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Introduction. In the development of sensor and measuring systems, modern small-sized household appliances, 
cellphones and wireless sensor systems, powerful sources of energy are not required for monitoring and diagnosing the 
technical condition of various objects, but mobility and nonvolatility of the above devices are mandatory. Energy 
storage devices with piezo-active elements that directly convert the energy of mechanical vibrations into electrical 
energy are widely used to power this kind of apparatus. In [1-3], energy storage devices using piezoelectric generators 
under the action of mechanical loads are studied. 

If the system is in an alternating magnetic field created by permanent magnets mounted on rotating parts of the 
machine, then the piezomagnetic layer is deformed along with the piezoelectric element. Due to this, an electric current 
is generated. In [7,8], Guo-Liang Yu and others discussed theoretical models of multilayer magnetoelectric composites 
for the magnetoelectric response at resonant frequencies corresponding to vibrations of the bending and extensional 
modes. A theoretical model of magnetic energy collectors using a functionally gradient composite cantilever was 
analyzed to improve the collecting ability and adjust the resonance frequency in [9]. In [10], Y.F. Zhanga and others 


studied bifurcations, periodic and chaotic dynamics of a tetrahedral composite multilayer piezoelectric rectangular plate 
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with simple supports. In [11], bending and free vibrations of a magnetoelectroelastic plate with surface effects were 
studied. 

In this paper, axisymmetric vibrations of the device are considered. The effect of the radius and thickness of 
the piezoelectric and piezomagnetic plates on the frequency characteristics of the device (eigenfrequencies of resonance 
and antiresonance) and on the efficiency of the conversion of vibratory energy into electrical energy, which is 
characterized by an electromechanical coupling factor, is studied. As a method for solving the problem, the finite 
element method implemented in the ANSYS package is selected. 


Z - axis of Symmetry 
a 
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Fig. 1. Scheme of half axial section of piezoelectric generator 


The energy storage device under consideration is an axisymmetric design, which consists of a metal disk 
(substrate) with two piezoelectric layers glued on it (Fig. 1). The top layer is piezoelectric; the bottom one is 
piezomagnetic. The flat surfaces of the piezoelectric layer are covered with electrodes that are connected to an external 
electrical circuit, or one of the electrodes is free while the other is set to zero electric potential. The piezomagnetic layer 
is affected by an alternating magnetic field according to the harmonic law; the outer radius of the substrate is pivotally 
fixed. The mathematical model of transverse steady-state vibrations of the described construction is the boundary-value 
problem of the linear theory of piezo-magnetoelectric elasticity [4]. The ANSYS package, which implements a 
piezoelectric model, is selected as a solution tool. 

The boundary value problem for a piezo-magnetoelectric body consists of a system of equations and boundary 
conditions [4]: 


V-o+pf=pu, V-D=0,, V-B=0 (1) 
o=c:e—-e -E-h’-H 
D=e:6e+K-E+a-H (2) 


B=h:e+o' -E+p-H 
e=—(Vu+(Vu)'), E=-Vo, B=-V6o, (3) 
where o and é€ are tensors of mechanical stress and strain; D and E are vectors of electric induction and electric field 
strength; B and H are vectors of magnetic induction and magnetic field strength; p is density of the material; c is 
tensor of elastic modules; e is tensor of piezoelectric modules; h is tensor of piezomagnetic modules; « is tensor of 
dielectric permeabilities; a is tensor of magnetoelectric modules; uw is permeability tensor; f is vector of the density 
of mass forces; 6, is volume density of electric charges; u is displacement vector; ~ and 9 are electric and magnetic 
potentials. 
Boundary conditions are specified for the mechanical, electrical, and magnetic components of the fields. 
Mechanical boundary conditions. Let the surface S consist of two parts T, andI’,, so thatS =T, UTl,, 
moreover, AT, =2. 
u=U onTI,, n-o=p onl,. (4) 
Electrical boundary conditions. Let the surface S consist of two parts T, andI,, so thatS=T,UT,, 
moreoverl, 01, =. 
p= (x,t) on T;, n-D=—o, on I,, (5) 
where o, is surface-charge density. In addition, if the electrodes are connected to an external circuit, two conditions 
should be added: 
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Ol, =v, [[n-DdS=I, (6) 


Sr 
where s, is electrode area; v is unknown potential which is found from the second condition; / is current. 
Magnetic boundary conditions. Let the surface S consist of two parts I; andI,, so thatS =I; UD,, moreover 
rPeor., =O. 
d=(x,t) on T,;, n-B=o, on I, (7) 
where o, is density of free surface currents along the boundary. 
For the elastic layer, the first equations in the system (1)-(3) are used; the components of the displacement 


vector u are unknown. For the piezoelectric layer, the second equations are added to them; u and the electric potential 


© are unknown. For the piezomagnetic layer, the third equations are added to the first equations, and the magnetic 
potential is added to the unknown displacements. In this case, the relations (2) are transformed through zeroing the 


corresponding constants. 

Finite element modeling. The computer model of the device is built in the ANSYS finite element package. 
The metal substrate (steel) has thickness fA and radius r. The piezo-active layers consist of one piezoceramic and one 
piezomagnetic plates which are polarized in thickness, have thickness h, and radius r, (Fig. 1). 

In the finite element model, devices for the metal and piezoelectric layers are used as finite elements PLANE42 
and PLANE13, respectively. In this work, the piezomagnetic layer is modeled by the finite element PLANE13, in which 
the piezoelectric properties of the material are replaced by piezomagnetic ones. This can be done for two reasons: the 
piezoelectric layers do not contact each other; qualitatively, the equations for the electric and magnetic potentials 
coincide. 

The properties of the piezoelectric layers used in the calculations are presented in Tables 1—2: piezoceramics 
— PZT-4, piezomagnetic material — CoFe204 [5, 6], adhesive layers are not taken into account. 

Table 1 


Material properties of PZT-4 piezoceramics 


e e e 
C",Hpa | C*%,Hpa | C..,Hpa | Ci, Hpa an al i | ome: | aw a 


77.8 74.3 iis 25.6 730 635 


€, = 8.85x10 F/m, PZT-4 piezoceramic density p = 7500 kg/m’ 





Table 2 
Material properties of CoFe,O, piezomagnetic element 


Ci. Cae Gas oe Cra Q;,, N/A Q;;, N/A O;;, N/A Mir» M3 > 


H Hpa N s/KI’ Ns/KI 


Hpa Hpa Hp 
2 


a pa m m m 
286 269.5 | 45.3 | 580.3 699.7 5.9x10" a0 





CoFe,O, piezomagnetic element density is p = 5290 kg/m° 
The elastic properties of the isotropic substrate material are characterized by the Young's modulus F and the 


Poisson's ratio v, E = 200 Hpa, v = 0.29 , density p = 7860 kg/m’ (steel) were used in the calculations. 


To achieve high accuracy in the calculations, the size of the final element of the metal layer is set to the value 


not higher than 1/5 its thickness; the size of the final element of the piezoelectric layers is automatically set. The finite 


element grid of the energy storage device 1s shown in Fig. 2. 
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Fig. 2. The finite element grid of the energy storage device. 
Surfaces 1, 4 are outer surfaces of piezoelectric layers; 
surfaces 2, 3 are contact surfaces between piezoelectric layers and metal layer. 


When analyzing the natural vibrations of a piezo-electromagnetic energy storage device, it is assumed that the 
following mechanical and electromagnetic boundary conditions are met. At the left end, symmetry conditions are 
specified, the right end is pivotally fixed (Fig. 1). To calculate the resonance frequency of the device, electric potentials 
are set on surfaces 1, 2; on surface 3, the magnetic potential is set; the magnetic-flux density is set on surface 4. In the 
case of calculating the antiresonance frequency of the device, the same boundary conditions as for calculating the 
resonance frequency are set on surfaces 1-4. However, the electric potential on surface 1 is unknown and is found from 
the condition (6). 

Numerical results. The natural vibrations of a piezoelectric element are considered, the radius of which is 
r,—9.8 mm, the thickness of the piezoelectric layers is h,=0.5 mm, the radius of the substrate is r=10 mm, and the 
thickness of the substrate is h=0.1 mm. Table 3 presents the first three natural frequencies of antiresonance. 


Table 3 
Antiresonance natural frequencies 


Antiresonance 


Distribution of vertical displacement on deformed state and 
distribution of electric and magnetic potential on undeformed sample 


frequency 
kHz 
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Antiresonance 


Distribution of vertical displacement on deformed state and 
distribution of electric and magnetic potential on undeformed sample 


frequency 
kHz 


115,372 





Table 3 presents the distribution of vertical displacement, electric and magnetic potential while the latter are 
presented in a single scale, here, the display of the magnetic potential is not visual. To build an approximate theory of 
calculating the vibrations of the considered bimorph, an analysis of the stress-strain state, electric and magnetic fields, 
which shows that for lower bending modes, there can be accepted hypotheses about their distribution corresponding to 
the bending of the plates. 

In the results of numerical calculations presented below, the dependence of the eigenfrequencies of resonance 
and antiresonance, and the electromechanical coupling factor on geometric parameters is investigated. 

When the number of alternating piezoelectric and piezomagnetic layers becomes sufficiently large, it 1s 
possible to use an approach based on the effective properties of the piezo-magnetoelectric composite [5, 6]. In this case, 
all equations of the problem (1) - (3) are used. 

The thickness value of the piezoelectric layers h, varies in the range of 0.3 + 0.7 mm, and the radius value r, — 


in the range of 6.8 + 9.8 mm. 


fal 





Fig. 3. Natural frequencies: a) resonance frequency, 
b) antiresonance frequency 


Fig. 3 shows the natural frequency dependences on the thickness h, and the radius r, of the piezoelectric 
layers. Fig. 3 shows that the values of the natural frequencies increase with increasing the radius. 

Fig. 4 shows the dependence of the electromechanical coupling factor on the thickness h, and the radius r, of 
the piezoelectric layers. Fig. 4 shows that the value of the natural frequency increases with increasing the radius of the 
piezoelectric layers r,, but decreases with increasing thickness h, of the piezoelectric layers. 
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Fig. 4. Electromechanical coupling factor 


Conclusion. An axisymmetric finite element model of an energy storage device based on round plates in the 
ANSYS package is considered. The active elements of the energy storage device are the piezoelectric and 
piezomagnetic plates fixed on a metal plate. The effect of the geometric characteristics of the piezoelectric layers under 
certain restrictions on the device dimensions, on the electromechanical coupling factor, which shows the efficiency of 
the energy storage device, is studied. The calculation results presented in the paper allow selecting rational sizes of 
piezoelectric elements operating at a certain frequency and having the greatest efficiency. This analysis shows that the 
maximum value of the electromechanical coupling factor is achieved when the thickness and radius of the piezoelectric 
layers take the largest and smallest values, respectively, within the considered limits. 
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